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The structures of the ground and excited states of 12 Be were studied with antisymmetrized 
molecular dynamics. The ground state was found to be a state with a developed 2a core with 
two neutrons occupying the intruder orbits. The energy levels of the newly measured spin-assigned 
states were described well, except for the lj~ state. The calculations indicated that many exotic 
cluster structures appear in the low-energy region. The widths concerning a and 6 He decays were 
discussed by using reduced width amplitudes. 
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Owing to the progress of experimental techniques, information concerning the excited states of light unstable nuclei 
has rapidly increased. Recently, exotic clustering in the light unstable nuclei has become one of the attractive subjects 
, in 3experimental and theoretical research. Since, in light stable nuclei, it has already been known that clustering is 
one of the essential features of nuclear dynamics, not only in excited states, but also in ground states, it is natural 
Q\ \ to expect cluster features in light unstable nuclei. Pioneering theoretical studies have suggested the development of 
cluster structures with a 2a core in Be and B isotopes 0,|,|,||J|,@,|,Ul^Jl^,|l|l . Especially, highly excited states 
with developed cluster structures in 10 Be and 11 Be have been studied by microscopic calculations Hfl^l 10 Xlf| . 
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In the case of 12 Be, the existence of cluster states was suggested in experimental measurements of the excited 
states |^3| . In recent experiments of 6 He+ 6 He and 8 He+ 4 He breakup reactions |p^ , p^| , many new excited states were 
discovered above the threshold energies. Some of the states are candidates of exotic cluster states, since the measured 
— . ■ spin-parities of those excited states indicate a rotational band with a large moment of inertia. It is an interesting 
subject to investigate clustering aspects in 12 Be. Although the molecular states in 12 Be were theoretically suggested 
with a potential model with He clusters [jl6|| , they have not yet been studied by microscopic calculations. In the recent 
investigations with the Generator Coordinate Method of a coupled-channel two-body cluster model, using microscopic 
6 He+ 6 He and 8 He+ 4 He wave functions, a strong mixing of both configurations was predicted Jl7| . 

The ground and the low-lying states of 12 Be present us with other attractive subjects concerning the vanishing of 
a neutron magic number, 8. The vanishing of the neutron magic number is already known in a neighboring nucleus, 
11 Be. The vanishing in 12 Be was predicted in the early theoretical works |l8[. The abnormal configurations((s<i) 2 ) in 
the low-lying states were experimentally studied by 10 Be(i,p) 12 Be reaction Jl9[. The abnormal structure of the ground 
state was also suggested in the analyses of the beta-decay strength from the ground state of 12 Be into 12 B |TT-j.[TTi|.pl~i| . 
Moreover, the vanishing in 12 Be was supported by a recent measurement of the spin-parity for a low-lying 1 — state 
plf . A new low-lying 0^" state H|j is also important to solve the inversion mechanism. In the recent theoretical 
studies, the structure of low-lying states were microscopically described by a molecular-orbital model by N. Itagaki 
et al. Q9|, and by a method of antisymmetrized molecular dynamics by one of the authors (Y.K .) |12[ |. They were also 
discussed by F. M. Nunes et al. within a three-body (n+n+ 10 Bc) model with core excitation [p3[. 

Our aim is to conduct systematic research of the structures of the ground and excited states of 12 Be based on 
microscopic theoretical calculations while focusing on clustering aspects. First of all, we have studied the systematics 
of the level structures, including the experimentally observed excited states. We have searched for cluster and non- 
cluster states to solve the following problems concerning clustering aspects in 12 Be. Do their cluster structures appear 
in Be ? If they do appear, what are the characteristics of their structures in unstable nuclei compared with those 
of stable nuclei ? The roles of the valence nucleons in the cluster states are interesting problems. We have also 
investigated the mechanism of the development and breaking of clustering. 

The important point is that the theoretical approach should be free from model assumptions, such as the stability 
of the mean-field and the existence of inert cores or clusters, because we have to describe various structures covering 
developed cluster structures as well as shell-model-like structures in the ground and excited states. It is difficult 
to study developed clustering in excited states with such mean-field approaches as the traditional shell model and 
the Hartree-Fock model. The cluster structures of the excited states of 9 Be and 10 Be were successfully explained by 
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cluster models pPM Eg , M while assuming a 2a core and surrounding neutrons. However, we think that a description 
based on cluster models simplifying the system as a 2a core with valence neutrons or 2-He clusters is not sufficient 
for a systematic investigation of 12 Be because of many valence neutrons in the system. We have applied a theoretical 
approach of antisymmetrized molecular dynamics (AMD). The AMD method has already proved to be a useful 
theoretical approach for the nuclear structure P,|5|,[7| pl]p5|j2^ | . Within the AMD framework, we do not need such 
model assumptions as inert cores, clusters, or axial symmetries, because the wave function of the nuclear system is 
written by Slater determinants, where the spatial part of each single-particle wave function for a nucleon is expressed 
by a localized Gaussian wave packet. Due to the flexibility of the AMD wave function, we successfully described 
the structure changes between shell-model-like states and cluster states as functions of the neutron number in light 
unstable nuclei. Owing to progress in computational power, it has become possible to study excited states by extended 
AMD calculations. This method is based on variational calculation after spin-parity projection (VAP calculation) 
within the framework of AMD, which has already been confirmed to be powerful for studying the excited states of 
light nuclei, as shown in studies of the stable nucleus 12 C pq |, and also unstable nuclei( 10 Be and n Be) ]^,pT|. The 
authors and their collaborators succeeded to describe various structures of the excited states and to reproduce many 
kinds of experimental data for nuclear structures of these nuclei with the AMD method. 

In the present work, the structures of the ground and excited states of 12 Be were analyzed by performing a variational 
calculation after spin-parity projection based on the AMD method. In the next section (Sec. ||), we explain the 
formulation of AMD for a study of the nuclear structure of excited states. The adopted effective interactions are 
explained in Sec. III. In Sec. [TV|, we present the calculated results concerning such observables as the energy levels, 



radii and (3 decays as well as the El and E2 transitions compared with the experimental data. In the discussion(Sec. 
|v|), the intrinsic structures and the rotational band structures are described. The single-particle behavior of the 
valence neutrons is analyzed. We discuss the systematics of the development of cluster states in Be isotopes, which 
can be classified according to the neutron orbits surrounding the 2a core. We investigate the inter-cluster motions 
between He clusters extracted from the obtained 12 Be wave functions, and calculate the decay widths for the 4 He and 
6 He channels by using the method of reduced width amplitudes. 



II. FORMULATION 

The formulation of AMD for a nuclear structure study of ground and excited states is explained in |^,^,^6|. In 
particular, the formulation of the present calculations is same as that described in Ref. 0. In this section, I briefly 
review the formulation. 

The wave function of a system is written by a superposition of the AMD wave functions^^/!))- The AMD wave 
function of a nucleus with a mass number A is a Slater determinant of Gaussian wave packets. The ith single-particle 
wave function is a product of the spatial wave function, the intrinsic spin function and the iso-spin function. The 
spatial part is presented by variational complex parameters, Xii, X2i, X3i, which indicate the center of the Gauusian 
wave packets. The orientation of the intrinsic spin part is expressed by a variational complex parameter ^, and 
the iso-spin function is fixed to be up(proton) or down(neutron) in the present calculations. Thus, an AMD wave 
function is expressed by a set of variational parameters, Z = {X„i, £j} (n = 1,2,3 and i — 1, • • • , A), which indicate 
the centers of Gaussians of the spatial part and the spin orientations of the intrinsic spin part of the single-particle 
wave functions. 

When we consider a parity-eigen state projected from an AMD wave function, the total wave function is written 
by two Slater determinants. In the case of a total-angular-momentum eigen state, the wave function of the state is 
represented by the integral of the rotated AMD wave functions. The expectation values of a given tensor operator for 
the total-angular-momcntum projected states are calculated by evaluating the integral with a sum over mesh points 
of the Euler angles. 

In principle, the total wave function can be a superposition of independent AMD wave functions. We can make the 
superposition of the spin-parity projected AMD wave functions (P^^^amd) as follows: 

$ = cPl± K ,<$> AMD (Z) + c'Pi± K ,$ A MD{Z') + ■■■. (1) 

We performed a variational calculation for a trial wave function to find the state which minimizes the energy of the 
system, 

(mm (9] 

by the method of frictional cooling, which is one of the imaginary time methods. Regarding the frictional cooling 
method in AMD, the reader is referred to two papers ||J25|]. 
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In order to obtain the wave function for the lowest J+ state, we varied the parameters Xj and £»(i = 1 ~ A) to 
minimize the energy expectation value of a spin-parity projected AMD wave function, $ = P M %t§ am x>(Z) , by using 
the frictional cooling method. That is to say, we performed energy variation after spin-parity projection(VAP) for an 
AMD wave function. 

With the VAP calculation for the J+ eigen state, ^^(Z) = P^-^amdC^), with an appropriate K', we obtained 
a set of parameters, Z = Z^, which presents the wave function of the first J+ state. In order to search for the 
parameters Z = Z^+ of the nth J+ state, the wave functions were superposed so as to be orthogonal to the lower 
states. The parameters Z„ for the nth J state were provided by varying Z so as to minimize the energy of the 
wave function orthogonalized to the lower states. 

After the VAP calculation of the states for various J, n and ir = ±, we obtained the optimum intrinsic states, 
$y4M_D(Z^), which approximately describe the corresponding J£ states. In order to obtain more precise wave func- 
tions, we superposed the spin-parity eigen wave functions projected from all of the obtained intrinsic states. Namely, 
we determined the final wave functions for the J+ states by simultaneously diagonalizing the Hamiltonian matrix, 

{PF K ,*AMD(Z*~ i )\H\P&%,*AMD(Z%* i )), and the norm matrix, (p£± ,^ AMD {Z J n f *)\Pj i ± IC „*AMD(Z$ r 3 )), with 
regard to (i, j) for all of the obtained intrinsic states, and to (K 1 , K"). Compared with the experimental data, such as 
the energy levels and E2 transitions, the expectation values were calculated with the final states after diagonalization. 

III. INTERACTIONS 

The adopted interaction is the sum of the central force, the spin-orbit force and the Coulomb force. The central 
force is chosen to be the MV1 force of case 3 [P7| , which contains a zero-range three-body force, V^ 3 \ as a density- 
dependent term in addition to the two-body interaction, V^ 2 \ of the modified Volkov No.l force. The spin-orbit force 
of the G3RS force || is adopted. 

IV. RESULTS 

The structures of the excited states of 12 Be were studied based on the VAP calculation within the framework of 
AMD. In this section, we present theoretical results, such as the energy levels, radii, and transitions (/3, E2, El and 
EO), while comparing them with the experimental data. A detailed analysis of the structures of the states is given in 
the next section. 

The adopted interaction parameters in the present work are those used in Refs. f^Jll[], which reproduced the 
abnormal spin-parity 1/2+ of the ground state of 11 Be. Namely, the Majorana, Bartlett and Heisenberg parameters 
in the central force are m — 0.65, b = h = 0, and the strength of the spin-orbit force is chosen to be ui = —uu = 3700 
MeV. The width parameter (v) was chosen to be 0.17 fm~ 2 , which gave a minimum energy of 12 Be in a simple AMD 
calculation without the spin projection. 

The wave functions for the lowest J+ states were obtained by the VAP calculation of P'1^ k ,^amd by choosing 
(J ± ,^')=_(0+,0), (2+,0), (4+,0), (6+0), (8+0), (1",1), (2" 1), (3~, 1), (4~ 1), (5~ 1) (6",1), (3+,2), (5+2), 
(7 + ,2), (0 ,0), (1+, 1), (l + ,0). After obtaining the lowest states(J{ t ), we calculated the second and third J ± states 
(0+, + , 2+, 2+, 4+, 6+, LJ ) with the VAP calculation of the superposed wave functions orthogonal to the obtained 
lower J states. The obtained AMD wave functions are considered to approximately describe the intrinsic states of 
the corresponding J+ states. The final wave functions of the J+- states were determined by superposing the spin-parity 
eigen states projected from these obtained AMD wave functions so as to simultaneously diagonalize the Hamiltonian 
matrix and the norm matrix. In principle, the number of the superposed AMD wave functions is 24, which is the 
number of calculated levels. In the present calculations, we omitted several of them which are not important, just 
to save computational time. Namely, we diagonalize the positive-parity(negative-parity) states projected from 22(18) 
AMD wave functions. 



A. Energies 

The theoretical binding energy of 12 Be was found to be 61.9 MeV(60.4 MeV) after(before) diagonalization, which 
underestimated the experimental value, 68.65 MeV. We found that the binding could be improved by changing the 
Majorana parameter (rn) of the interaction to a smaller value. For example, a theoretical value of 66.1 MeV was 
obtained by using m — 0.62 with a simple VAP calculation before diagonalization. Nevertheless, we adopted the 
parameter m — 0.65 in the present work, because this value reproduces the parity inversion of n Be 0, which should 
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be important for describing the ground-state properties of the neighboring nucleus, 12 Be. We checked that the change 
in the m parameter has no significant effect on the excitation energies, £?(0f ) and E(8f), of 12 Be, at least. In fact, 
the calculated excitation energies of S(0f ) before diagonalization are 3.0 MeV and 2.8 MeV, in the cases of m = 0.62 
and m = 0.65, respectively. 

The theoretical level scheme is shown in Fig.0. The calculations suggest that many excited states appear in the 
low-energy region. This is the first theoretical work which systematically reproduces the experimental energy levels of 
all the spin-assigned states, except for the 1~ state. By analyzing the intrinsic AMD wave functions of the states, we 
can consider that there exist rotational bands (K n =0f , Of, Of and which consist of the following states: {Of, 
2+, 4+ , 6f and 8f }, {Of , 2+}, {Of , 2+, 4+ and 6f } and {lj~, 2^, 3^, 4f, 5f}, respectively 

The energy-spin systematics for positive-parity states with natural spins are shown in Fig. |^. It is surprising that 
the newly observed levels, 4 + at 13.2 MeV and 6 + at 16.1 MeV [Q, correspond well to the 4 f and 6 J states obtained 
in the present results. An interesting point is that these excited states belong not to the yrast band, = Of, but 
to the new excited = Of band with developed cluster structure. The theoretical results predict the existence of 
many positive-parity states belonging to the K 7 ' = Of and = Of bands in the lower energy region. Even though 
12 Be has a neutron magic number of 8, the calculated ground K v = + band has a large moment of inertia, because 
the band head Of state is not the ordinary state with a closed neutron p-shell, but an intruder state with a developed 
cluster structure. It has a prolate deformed structure, which is dominated by 2p-2h configurations, and reaches the 
band terminal at the 8f state, because the J = 8f state is the highest-spin state in the 2hu> configurations. The 
third 0+ band (K * = 0J ) was found to be represented by other 2p-2h configurations than those in the K v = Of band. 
The K v = 0^ band has an extremely large moment of inertia because of a remarkably developed 6 He+ 6 He clustering, 
and reaches the band terminal 6f state accompanying the spin-alignment of nucleons in the high-spin region. On 
the other hand, the main components of the Of and 2f states are of the OTiw configurations, with the closed neutron 
p-shell, and constitute the = Of band. The band head Of state of this band, which was theoretically predicted to 
appear just above the intruder ground band by Itagaki et al. B and Kanada-En'yo et al. |32j, was recently discovered 
in an observation of coincidence gamma rays by Shimoura et al. J23|. 

The theoretical excitation energy of the lj~ state is larger than the measured 1~ state at 2.68 MeV Ee. In spite 
of the overestimation of the excitation energy, we regard this lj" state in the K* = 1~ as the l - state at 2.68 MeV 
because of the large El transition strength, as shown later. In order to improve the excitation energy of this 1~ 
state, it is important to take the mixing of the other configuration into account. For example, by mixing a 1 _ state 
with K v = _ , which is obtained at slightly higher energy than the 1~ state with K n = 1~ in VAP calculations, the 
1^ state gains about 1 MeV. Another reason for the overestimation may be connected to the fact that the present 
interactions give a too small value for the effective energy difference between the S\/2 and d§/2 orbits, which was found 
in the calculations of 11 Be in Ref. JO]]. 

B. Radii 

The theoretical values of the root-mean-square radii of the density distributions of point-like nucleons, protons, 
and neutrons are listed in Table | together with the experimental matter radius j2j| deduced from the reaction cross 
sections. Because of deformations, the proton and neutron radii in the ground state are larger compared with those 
in the Of state. Since the theoretical values concerning the difference Ar = r n — r p between the proton and neutron 
radii are about 0.3 fm, 12 Be is a candidate of the neutron skin nucleus. Compared with the experimental data, the 
present result for Of is larger because the present parameter m = 0.65 is considered to be too large to quantitatively 
reproduce the radii of p-shell nuclei. 

The proton and neutron densities as functions of the radius are shown in Fig. |^. In both the Of and Of states, the 
excess neutrons enhance the neutron density in the surface region. In the Of state, due to cluster development, the 
proton density increases in the surface region while it decreases at the center of the nucleus. It should be noted that, 
even if a nucleus has a neutron halo (the long tail of the neutron density in the outer region), the details of the halo 
structure are not expressed in the present model space because of a limitation of the Gaussian forms in AMD wave 
functions. 



C. (3 decay strength 

The strength of the f3 decay from 12 Be(0f ) to 12 B(1 + ) provides helpful experimental evidence for breaking of the 
neutron p-shell closure in 12 Be(0f ). T. Suzuki et al. suggested that weak (3 decay is experimental evidence for an 
admixture of non-0?ia; configurations in 12 Be(0f ), because the f3 decays from a normal p-closed state of 12 Be(0f ) must 
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FIG. 1. Excitation energies of the levels in 12 Be. The experimental data are from the Table of Isotopes and Refs. 

iliilBJiiEl. 
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FIG. 2. Energy-spin systematics for positive-parity states in Be. The excitation energies of the natural spin states are 
plotted as functions of J(J + 1). The circle symbols are the theoretical results, while the cross symbols are the experimental 
data of the spin-parity assigned states, which are taken from the Table of Isotopes and Refs. Jl4| , p2"[ . 
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TABLE I. The root-mean-square radii of the density distribution of point-like nucleons, protons, and neutrons. The experi- 
mental matter radius deduced from the reaction cross sections is taken from Ref. ||2Sf| . 





matter 


proton 


neutron 


12 Be(0+) 


2.85 fm 


2.67 fm 


2.94 fm 


12 Be(0+) 


2.75 fm 


2.56 fm 


2.84 fm 


exp. 


2.59 ± 0.06 
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FIG. 3. Densities as functions of the radial coordinate (r). The densities of the protons and neutrons in the 0^ state (solid) 
and in the 0^" state(dashed) are shown. 
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TABLE II. The strength of f3 decays. B(GT) is defined as (<tt)| 2 . The experimental data are taken from pi 
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final 
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12 BefO + 0) 


12 B(1+ 0) 


0.56 
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before diagonalization 


after diagonalization 


ii! Be(0+) 


-B(l+) 


0.2 


0.9 


12 Be(0+) 


12 B(1+) 


3.0 


2.1 


12 Be(0+) 


12 B(1+) 


0.05 


0.07 


12 Be(2+) 


12 B(1+) 


0.004 


0.1 


12 Be(2+) 


12 B(1+) 


0.5 


0.4 



be stronger. We calculated the Gamow- Teller transition strength, B(GT) = |(crr)| 2 , where the wave function for the 
daughter state, 12 B(1+), was obtained by a VAP calculation with ( J ± , K') = (1+, -1). The B{GT) values are given 
shown in Table [fl| The present result, B(GT) — 0.9, is almost as small as the experimental data, B(GT) = 0.59, 
because the component of the 2hu> configurations in the parent 12 Be(0^~) makes the transition matrix element of the 
Gamow- Teller operator to be small. In fact, the original Of state of 12 Be obtained by a simple VAP calculation before 
diagonalization has very weak (3 transitions as B(GT : 12 Be(0f ) — ► 12 B(lf)) = 0.2. On the other hand, the decay from 
Of is strong as B(GT : 12 Be(0^) — > 12 B(lf)) = 3.0 before diagonalization, because the parent state has the ordinary 
Ofiui configuration. After state mixing with diagonalization, the GT strength from 12 Be(0^~) is distributed to that 
from 12 Be(0^~). As a result, after the diagonalization the B(GT : 12 Be(0^~) -^ 12 B(lf)) value increases to 0.9, while 
B(GT : 12 Be(0^) — > 12 B(1^)) decreases to 2.1. According to the present results, the GT transitions from 12 Be(0^) are 
weak, while the B(GT : 12 Be(2+) ^ 12 B(1+)) is largest among the 2 + states because the 2% state is the other Ohui 
state belonging to the K 71 — Oj band. 

The present results concerning the weak GT decay from the ground state is consistent with the discussion given in 
Refs. 020). 



D. E2,E1,E0 transition strength 



Although the data concerning the E2 transition strength provide good information about proton deformations, we 
must take care of the following points in the analysis of light neutron-rich nuclei. Firstly, it is dangerous to assume the 
same deformation of the proton density as that of the neutron density in light unstable nuclei. Secondly, the classical 
relations between the intrinsic deformation (Qo) and the observables, which are often used in a simple analysis of 
heavy nuclei, do not necessarily work in light nuclei. Therefore, it is necessary to analyze the transition strength based 
on a microscopic calculation. As shown in previous AMD studies [p|,p6|j7|l, the experimental Q-moments and B(E2) 
values of light nuclei were reproduced well by using bare charges, because of the advantage of the AMD method, 
which can directly express the proton and neutron deformations. 

Table III gives the theoretical B(E2), B(E1) and B(E0) values calculated by the VAP calculation after diagonal- 
ization. The intra-band E2 transitions in the = Of band are strong, as seen in B{E2\2\ — > 0^) = 14 e 2 fm 4 , due 
to the deformed intrinsic state. On the other hand, B(E2) is smaller in the transition between the 2^ and Oj" states 
in the = Oj" band, which have rather spherical shapes compared to those in the X* = Of band. Although the 
deformations of the states in the K* = 0% band are extremely large in the present results, B(E2; 4^ — > 2^) and 
B(E2; 2\ — ► O3 ) are not as large as those in the K v = Of band, because the intrinsic structure varies with increasing 
the total spin (J) along the 0^~ band. The B(E2) values regarding the transitions between different bands are smaller 
compared with those for the intra-band transitions. 

Compared with the experimental data, B(E2; 2^ — > 0f) = 10.5±1.1 (e 2 fm 4 ) of 10 Be, the present results predict a 
larger B(E2) of 12 Be as B(E2; 2\ — > Of ) =14 e 2 fm 4 . We should point out that the E2 transition strength is sensitive 
to the mixing ratio with the 0j~ and 2\ states as well as deformations. Since the 0^" and 2^" states are suggested to 
exist just above the Of and 2f states, we should more carefully investigate the energy difference and the state mixing 
between the = Of and — 0^ bands before making conclusion. 

The El transition strength was recently measured as B(El;0f — > 17) = 0.05 e 2 fm 2 pi] , which is rather large 
compared with other light nuclei. The calculated result, B(El;0f — > 1^)=0.02 e 2 fm 2 , reasonably agrees with this 
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TABLE III. Theoretical results of E2, El and EO transition strength, B(EX). B(E0) is defined as | < /|i+pr 2 |i > 



transitions 


band(_ftT' r ) 


Mult. 


present 


i2 Be;2+ -> 0+ 


in 0+ 


E2 


14 (e a fm 4 ) 
8 (e 2 fm 4 ) 
8 (e 2 fm 4 ) 


12 Be;2+ -> 0+ 


m 0+ 


E2 


12 Be;2+ 0+ 


in0 3 + 


E2 


12 Be;4+ -> 2+ 


in 0+ 


E2 


14 (e 2 fm 4 ) 
10 (e 2 fm 4 ) 
5 (e 2 fm 4 ) 

5 (e 2 fm 4 ) 

6 (e 2 fm 4 ) 


12 Be;4+ -> 2+ 


in 0+ 


E2 


12 Be;4+ - 2+ 




E2 


12 Be;4+ -> 2+ 




E2 


12 Be;4+ -> 2+ 




E2 


i2 Be;0+ -» 17 




El 


2 x 10~ 2 (e 2 fm 2 ) 
2 x 10 -2 (e 2 fm 2 ) 


12 Be;0+ -> 17/ 




El 


12 Be;0+ 0+ 




EO 


1.7 (e fm 2 ) 



experimental data. Before diagonalization, the deformed 0+ state possesses the strength of the El transition from 
the 17/ state, while the transition to the spherical 0+ state is weak. After diagonalization, the strength in the 0+ state 
distributes in the 0+ state due to state mixing between the 0+ and the 0+ states. In other words, the rather large 
B(E1; 0+ — * 17/) is caused by a deformati on of the ground state. If the mixing of the K v = CP component in the 
state is taken into account as described in IV A, theoretical value of the El strength become large as B(E1; 0+ — * 17/) 
=0.14 e 2 fm 2 , because the K 1 * = 0~ component has a similar intrinsic structure as that of the 0+ state. 



V. DISCUSSIONS 



In the present calculations, many deformed rotational bands appeared in 12 Be. It was found that 2a core is formed 
in most of the states as well as in other neutron-rich Be isotopes: 10 Be and 11 Be [^Mp] . In this section, we consider 
the intrinsic structures of 12 Be while focusing on the clustering aspects and the behavior of the valence neutrons. The 
features of the intrinsic structures, the roles of valence neutrons and the inter-cluster motions in the cluster states are 
discussed. 



A. Systematics of rotational bands 



By analyzing the structures of the intrinsic states, we can classify the ground and excited states into rotational bands. 
We superpose the wave functions, Pm K $amd, projected from all of the intrinsic wave functions, $ J 4Mr)(Z^ , .' 7ri ), 
obtained by VAP so as to diagonalize the Hamiltonian matrix. Although the final wave function $j± ( 12 Be) for the J+ 

state after diagonalization is the linear combination of spin-parity projected AMD wave functions ^M^iMD(Z^'"'), 
we consider the AMD wave function before diagonalization as being the intrinsic state of the corresponding J+ state, 
because P^$amd(Z^) is found to be the major component of the final result $j±( 12 Be) of the J+ state, except 

, the squared amplitudes, ($j± ( 12 Be)|P^^$AMT>(Z;[+-) 2 are listed. The 

amplitudes are larger than 0.6, except for the 2+ and 2+ states. In the negative-parity state of 1^, the component 
of the original state becomes less than 0.6 after diagonalization. Considering that this state might be unstable, we 
neglect the 1^~ level in the present results. In the case of the 2+ state, the main component of the 2+ state is not the 
projected AMD state, P 2 / +$aa / .d(Z 2 +), obtained by a VAP calculation for J+ = 2+, but is the P^ 2 ^ AM d{'^\ + ), 
which is projected from the intrinsic state of the 3+ state. This means that a new 2 + state appears in addition to three 
2+ states obtained by the VAP calculation, Z 2 +, Z 2+ and Z 2+ . The energy of this new 2+ state, P^^amdCZ^), 
is slightly lower than that of Pj 1 ^,$^mo(Z 3 + ). As a result, the wave functions, $amd(Zi + ) and 3?amt>(Z 3 +), are 
regarded to be the intrinsic states of the 2+ and 2+ states, respectively, although the amplitude of the original 
VAP state, P^ q ^amd{^1 + ), hi the 2+ state is as small as 0.4 after diagonalization because of mixing the 2+ state 
with the 2+ state. The reason why the intrinsic state of the 2+ state is not obtained in the VAP calculation with 
{J ± ,K I ) = (2+,0) is because the 2+ state belongs to the K 7 * = 2+ band. 

By analyzing the intrinsic states, we find the rotational bands = 0+, 0+, 0+, 2+, 17/, which consist of the states 
(0+2+ 4+ 6+8+), (0+2+), (0+,2 + ,4+6+) (2+ 3+ 5+ 7+ ), and (17/, 27", 37/, 47, 57), respectively. We analyze the 



for the 2+ and 2+ states. In Table IV 
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TABLE IV. The amplitudes of the corresponding VAP wave functions Pjjfa &AMD('Zn/ Ti ) in the final wave functions 
($ J ±( 12 Be)) obtained after the diagonalization. Total intrinsic spins (S 2 ,) and (S 2 ) for protons and neutrons in the spin-parity 

projected states P^^AMcfZif') before the diagonalization are also listed. 



Jn after diagonalization 


(J^M-K^n*) in VAP 


<$ /± ( 12 Be)|P^±$ AMD (zC)> 2 


<s*> 


(Si) 


or 


(0+,0, 1) 


0.8 


0.1 


1.0 


9+ 


, u, i ) 


U.o 


U.l 


1 1 
1.1 


*t 


(4 + ,0, 1) 


0.9 


0.1 


1.2 




(6+,0, 1) 


0.9 


0.0 


1.5 


«t 


(8+,0, 1) 


0.9 


0.1 


1.6 


o 2 + 


(0+,0,2) 


0.6 


0.2 


0.1 


2 2 + 


(2+,0,2) 


0.7 


0.9 


0.1 


o 3 + 


(0+,0,3) 


0.7 


0.0 


0.5 


2+ 


(2+,0,3) 


0.4 


0.0 


1.0 




(4+,0,2) 


0.6 


0.0 


1.0 




(6 + ,0,2) 


0.6 


0.0 


1.2 


ir 




0.9 


0.1 


0.8 


2r 


(2 _ , 1, 1) 


0.9 


0.1 


0.8 




(3", 1,1) 


0.9 


0.1 


1.0 


4r 


(4",1, 1) 


0.7 


0.1 


1.1 


sr 


(5", 1,1) 


0.9 


0.5 


1.2 


2 3 + 


(3+,2, 1) 


0.8 


0.0 


2.1 


3? 


(3 + ,2, 1) 


0.8 


0.0 


2.1 


si- 


(5 + ,2,l) 


0.8 


0.0 


2.0 


7i- 


(7 + ,2, 1) 


0.9 


0.0 


1.2 


or 


(0",0, 1) 


1.0 


0.1 


1.7 


ii- 


(1+, 1, 1) 


0.9 


0.0 


1.3 


$ 


(1 + ,0,1) 


1.0 


1.9 


0.3 




(6", 1,1) 


0.9 


0.0 


2.0 
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parity-eigen components in the single-particle wave functions of the intrinsic states, and find that each rotational 
band is dominated by either the OTiuj, lfuv or 2fuv excited configuration. The states in the three positive-parity bands 
(A 71 " — 0^ , 0^, 2^ ) have the dominant 2Huj excited configurations, while the states in the = 0^ band are dominated 
by the 2hu> configurations. The first point is that the states in the ground band are not the neutron shell-closed states, 
but are the prolately deformed intruder states with developed clustering, even though this 12 Be nucleus has a neutron 
magic number of N = 8. This is consistent with discussions of the vanishing of the neutron magic number in 12 Be 
given by Refs. ]l8| , |T9| ,p|, pTpO| ] . The A x = Of band starting from the state terminates with the 8^~ state, because 
the J n = 8 + is the highest spin in the 2hui configurations. On the other hand, the present results predict that the 0^ 
and 2^ states with the neutron p-shell closed structures constitute the excited band, = 0^", just above the ground 
band. The second interesting point is that the newly observed levels, 4 + and 6 + states |TJ], are considered to belong 
to the higher band of the two bands(A' 7r = 0^ and K n = 0^ ), both of which are the 2Huj excited states with developed 
cluster structures. The K n = 2± band consists of the 2^,3^,5^,7^" states, which have the other 2hcu configurations 
. Although spin-parity eigen states, 4 + and 6 + , are able to be projected from the intrinsic states of this K n = 2+ 
band in principle, the 4 + and 6 + states in the if" = 2\ band can not be identified, because these states strongly mix 
with the 4+ , 4^, 6f and 6 J states in the = and A T — 0jj~ bands. We find a negative-parity band, A T = 1^, 
in which the states(l _ ,2 _ ,3 _ ,4 _ ,5~) are dominated by lhiu configurations. 



B. Intrinsic structures 



We discuss the features of the intrinsic structures, such as the deformations and clustering aspects. The density 
distributions of the intrinsic AMD wave functions are shown in Fig. ^. 

In most of the states, except for the 2 J, 5j~ and states, the densities of protons have dumbbell-like shapes, which 
indicate the formation of a 2a core. On the other hand, 2a core breaking occurs in the 2J, and states. For 
the quantitative discussions, we can estimate the degree of core breaking based on the non-zero values of the squared 
total- intrinsic spin of protons, (S 2 ), given in table [V. In most of the states, the values (S 2 ) almost equal to 0, which 
indicates only a slight breaking of the 2a core. However, as seen in the large (S 2 ) values in the 2\ and 5j~ states, the 
2a core structures somehow dissociate in these band terminal states of the K v = 0^ and l - bands because of the 
spin- alignment effect. On the other hand, the 1^ state has a quite different structure from that of the 2a-core state. 
It is a spin- aligned state with proton spin, S p ~ 1, where one of the 2-a clusters is completely broken. In fact, no 
dumbbell shape is seen in the proton density of this state (Fig. ^). The alignment of the intrinsic spins is the origin 
of the unnatural spin parity of this state. This \\ state is analogous to the 3/2~ state at 8.4 MeV in 11 Be, and also 
to the 1± state in 12 C, in which cluster breaking occurs due to the the aligned intrinsic spins. 

The deformation parameter of the proton density shown in Fig. |] is a useful quantity in order to quantitatively 
discuss the developments of clustering in the 2a core states. Here, we explain the cluster developments in relation to 
the neutron structure, such as the neutron deformation and neutron intrinsic spin. 

In the K 1 * = Of band, the states have a prolate deformation with the developed 2a core, as shown in Fig. ||. In the 
0^ , 2^ and 4^~ states, the 2a core well develops following prolate deformation of the neutron density, which is caused 
by the 2 neutrons in a longitudinal sd-orbit. The ground Of state has an approximately axial symmetric shape (Fig. 
^), which is considered to be caused by the symmetric neutron configuration described by a p 3 / 2 sub-shell closure 
and the axial symmetric sd-orbit with (X/i) = (20) symmetry in the SU3 limit. The non-zero total intrinsic spin of 
neutrons of the ground 0^ state shown in Table |W is considered to come from a sub-shell closure effect of the neutron 
p 3 / 2 -shell. As the total spin increases, the alignment of the intrinsic spins of the sd-shell neutrons further increases 
the total intrinsic spin of neutrons (Table |TV|). The aligned sd-shell neutrons make the neutron-density deformation 
of the nucleus to be smaller in the high-spin region, J > 6. As a result, the clustering become weak in the and 8f 
states, as found in the reduction of the proton deformation parameter (Fig. [|). 

In the A 7r = 0^" band with the dominant Ohiv configurations, the neutron deformation is smaller than that in the 
A 71 " = 0^ band. In the 0^ state, 4 neutron pairs make a small tetrahedron-shape structure, which is approximately 
equivalent to the p-shell closed structure. Because of the shell effect of the neutron p-shell closure, the clustering is 
smaller than that in the Of state. Although the intrinsic state of the band terminal 2\ state in this band has the 
prolate neutron density as seen in Fig. ||, the spin-parity J ± = 2 + eigen state projected from this intrinsic state 
is dominated by the neutron p-shell closure component. The total spin J = 2 in this 2\ state is composed of the 
aligned total-angular momentum of the p-shell protons. The overlap \{P^ K=0 ^AMD{^ { 2 + )\PlfK=o^AMD{'^2 + ))\ 2 is 
still large as about 0.6. The j-j coupling feature in this band terminal state is found in the non-zero value of the 
squared total-intrinsic spin of protons (Table [V) and also in the disappearance of the dumbbell- like shape of proton 
density (Fig. |). 

In the A 11 = 0^ band, the 0;j~ and 2\ states have extremely developed cluster structures, like 6 He+ 6 He. One of 
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FIG. 4. Density distribution of the intrinsic state $amd(Z^) before projection and diagonalization. The intrinsic system 
is projected onto a plane which contains the approximately longitudinal axis of the intrinsic state. The density is integrated 
along a transverse axis perpendicular to the plane. The densities for matter, protons and neutrons are presented at the left, 
middle and right, respectively. The deformation parameters /3 are written at the bottom of the figures. The expectation values 
of xx, yy, zz for the matter, proton and neutron densities are written below the figures. The size of the frame box is 10 fmx 10 
fm. 
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the reasons for the development of clustering is considered to be the orthogonal condition to the Of and 2f states in 
the lower band. The inter-cluster distance shrinks with increasing total spin, due to the spin alignment. In the Of, 
2f and 4if, low matter density regions between the clusters are seen (Fig. ||), while they disappear in the 6 J state. 
The intrinsic structure of the = Of band changes rather rapidly with an increase of the total spin, as can be seen 
in the decrease of the deformation parameter (3 (Fig. |J) and also in the spin alignment of the neutron intrinsic spins 
(Table IV). Eventhough the 6 He+ 6 He cluster structure changes due to the spin alignment in this K* = Of band, 
it is help to consider the spatial symmetry in the SU3 representation. Since the intrinsic state of the band head Of 
has (A/i)=(24) symmetry in the SU3 limit, which provides a K — band from J 77 = + — > 6 + , it is reasonable that 
K* = Of band terminates at the state. 

In spite of the mixing of the states and the spin alignments, all of the states in the K n = Of, Of and 2f bands 
are dominated by 2hu> configurations with 2 neutrons in the sd-shell. It is another interesting problem to search 
for the weak-coupling cluster structure, which is described by the relative motion between the 2- 6 He clusters. In 
pioneering work by M. Ito et al. [jLfj, weak coupling states with 2- 6 He clusters in 12 Be were studied. Although the 
calculation was not fully microscopic, because they did not make antisymmetrization of the neutrons between clusters, 
it is interesting that He-cluster states were suggested to appear near and above the threshold energies. The presently 
predicted Of state is a candidate of the weak coupling cluster state, because of the large relative distance between 
clusters. However, in the and 6^" states in the K 77 = Of bands, the components of the weak coupling cluster state 
may not be large, since the He-cluster structure becomes weak due to the spin alignments. Instead, it is expected 
that the molecular resonances with a weak-coupling cluster structure may exist above those spin-aligned states. We 
performed a simple VAP calculation of the higher states, and found that candidates of the weak-coupling states appear 
above those spin-aligned states with 2Tuo configurations obtained in the present calculations. We find a 4f state with 
a developed 6 He+ 6 He-cluster structure at a 0.6 MeV higher energy than the 4f state. Although the 6f state was 
unstable for cluster escaping in the present VAP calculations using m = 0.65, we obtained a 6f state with a developed 
6 He+ 6 He structure at a few MeV higher energy than the 6^" state by using m = 0.62. Such weakly bound states 
should be carefully investigated by taking care of the threshold energies and the stability for particle decays. 

The negative-parity states If, 2f , 3f , 4f and 5f in the K 77 = If band are deformed states with cluster structures. 
The deformed structures of the neutron density are made from the dominant IHlu excited configurations with one 
neutron in the sd-shell. The degree of prolate deformation of the neutron density in this band is smaller than the 
states in the K 77 = Of band. The smaller neutron deformation causes the weaker cluster development in the K 77 = If 
band compared with those in the = Of band. The If, 2f , 3f , 4f states have axial-asymmetric neutron structures, 
while the band terminal 5 _ state has approximately axial-symmetric neutron density. The reason for K 77 = 1~ of 
the If state can be naturally understood by a weak coupling picture of 11 Be core and a neutron. In the simple weak 
coupling picture, since the If state of 12 Be can be described by the coupling of a n Be(l/2~) core with (A/x) = (21) 
SU3 symmetry and a s^-orbit neutron, the K 77 = in the 12 Bc(l~) state come from the K 77 = 1~ in the n Be core. 
In the recent work with a three-body(n+n+ 10 Be) model |23|], in which the low-lying states of 12 Be were studied from 
the weak coupling picture, the excitation energy of the If state was predicted to be about 3 MeV. Another possible 
negative parity band is K 77 = 0~, which is a parity doublet band of a 8 He+ 4 He cluster structure, as is suggested 
in a He+He cluster model by P. Descouvement et al. [0. We found a l - state with the 8 He+a cluster structure 
by the VAP calculation for K 77 = 0~ at slightly higher energy than that for K 77 = 1~, Even though the mixing of 
this K 71 = 0~ component is important to lower the excitation energy of the If state as mentioned before in the 
major component of the If state is the K 77 — 1~ state. The improved results including the = 0~ configuration 
suggest the higher state(lf ) with the significant K 77 = 0~ component may appear at 9 MeV excitation energy, which 
is consistent with the prediction in Ref. | fj"7| . 

Next, we discuss the structures of the unnatural spin-parity states. As mentioned above, the unnatural spin-parity 
of the 1 f state is caused by the aligned intrinsic spins of protons. On the other hand, the unnatural spin-parity of the 
other states (3f , 5f , 7f, If and Of) originates from the aligned intrinsic spins of the s<i-shell neutrons. Alignments 
of the neutron intrinsic spins were found in the 2^~, 3^~ and 5f states, as seen in the expectation values of the squared 
total neutron spins, (S 2 ) ~ 2, presented in Table pVL The spin alignment is mainly caused by 2 sd-shell neutrons. 
The deformation and the density distribution of the intrinsic state of the 7f state are very similar to those of the 3f 
and 5f , though (S 2 ) in the 7f state is slightly smaller. The If and Of states with unnatural spin-parity have aligned 
intrinsic spins of neutrons, like the states in the K 77 — 2f band. Although a 6 He+ 6 He-like structure is found in the 
if state, one of the 6 He clusters is an excited 6 He with S n = 1 components. 

Here, we stress that it is essential to consider the orthogonality to the lower states as well as antisymmetrization 
of all the nucleons in a study of the excited states of 12 Be, because the present results suggest the existence of many 
2a core states in the low-energy region. Although these states are not necessarily written in terms of 6 He+ 6 He 
clusters, such lower states must have effects on the higher weak-coupling cluster states. The present results indicate 
the importance of treating degrees of freedom of the intrinsic spins for all the nucleons in a study of 12 Be. The reasons 
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FIG. 5. Sketches for the molecular orbits, (a) u-orbits, (b)-7r-orbits, and (c)5'-orbits surrounding 2a core. These molecular 
orbits are explained by linear combinations of the p-shell orbits around the a cores. 



can be summarized as follows. First of all, a closed p 3 / 2 -shell of neutrons plays an important role for the energy 
gain of the intruder ground state. Secondly, the alignments of intrinsic spins as S = 1 are necessary to describe the 
unnatural spin-parity states. In the third point, the spin alignment of the neutron intrinsic spins occurs in high spin 
states. 



C. Picture of molecular orbits 



As mentioned in the previous subsection, the 2a core structure appears in many of the states. The development of 
clustering is sensitive to the prolate deformation of the neutron density. We inspected the neutron structure based on 
the single-particle behavior of the valence neutrons surrounding the 2a core. We also considered the relation between 
the neutron orbits and the cluster development. We analyzed the single-particle wave functions of the valence neutrons 
in the intrinsic states. The detailed formulation of how to extract the single-particle energies and wave functions from 
an AMD wave function is described in Refs. Hf5J. 

According to the analyses by the single-particle orbits and the single-particle energies, the 2a core is composed 
of nucleons occupying the lowest 4 proton orbits and 4 neutron orbits. The higher 4 neutron orbits correspond to 
those of the valence neutrons surrounding the 2a core. By extracting the positive and negative components in the 
orbits, many of the single-particle orbits are found to be approximately parity-eigen states, as the amplitude of the 
dominant parity-eigen state is more than 70% in each orbit. Since the negative and positive-parity orbits of the valence 
neutrons are associated with the p-orbits and sd-orbits, respectively, we can classify the ground and excited states of 
12 Be in terms of the nfiuj excitation, where n is the number of the valence neutrons with the dominant positive-parity 
components. The K 7 * — 0^ band is approximately described by OHuj configurations, while the = 0^ , 0^~, 2+ bands 
are dominated by 2tuv configurations with 2 neutrons in sd-like orbits. On the other hand, the main components of 
the K" = li band are lhto configurations. 

The idea of molecular orbits surrounding a 2a core is helpful to understand the roles of the valence neutrons in 
neutron-rich Be isotopes. The molecular orbits in Be isotopes were suggested in a study of 9 Be with a 2a + n cluster 
model pi} ]. They assumed cr-orbits and 7r-orbits which are made from linear combinations of the p-orbits around the 
a cores (see Fig.g). This idea was applied to neutron-rich Be isotopes by Seya et al. a long time ago [|lj. In the 1990's 
von Oertzen et al. [3,31j revived this kind of research to understand the rotational bands of neutron-rich Be isotopes, 
and Itagaki et al. [£, ^described the structures of the low-lying states of 10 Be and 12 Be by assuming 2a core and 
valence neutrons in the molecular orbits. The formation of the 2a and valence neutron structures in neutron-rich Be 
isotopes was first guaranteed theoretically by the AMD calculation ||,||,^],f7||ll]] , where the existence of any clusters or 
molecular orbits was not assumed. In these AMD studies, the viewpoint of the molecular orbit was found to be useful 
to understand the cluster development in 10 Be and 11 Be. Therefore, it is an interesting problem whether the states 
of 12 Be can be described by the molecular orbits. 

In the present results for 12 Be, we find a new kind of molecular orbit besides the suggested 7r-orbit and cr-orbit. 
In the positive-parity orbits of the valence neutrons in 12 Be, two kinds of molecular orbits appear, both of which are 
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FIG. 6. Density distributions of the single-particle wave functions of the valence neutrons in the intrinsic wave functions of 
the Oj", O3 , 6^ and 3f states. The figures at the left (right) show the densities regarding the positive-parity components of the 
first(second) highest neutron orbits. The value P+ in each orbit indicates the squared amplitude of the contained positive-parity 
component. 
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associated with sd-orbits. The first one is the cr-orbit (Fig. ||a), while the second one is a quite new molecular orbit, 
shown in Fig. ||c. This orbit is the other positive-parity orbit made from a linear combination of the p-orbits around 
the a cores. As shown in Fig. [|c, the combined p-orbits in this orbit are perpendicular to those in the cr-orbit. We 
call this new positive-parity orbit a cT-orbit in the present paper, although it has (A/i) = (01) symmetry in the SU3 
limit, which is perpendicular to the so-called (5-orbit in the field of the molecular physics. In the case of 12 Be, the 
negative-parity orbit of the neutron surrounding 2a does not necessarily correspond to the pure molecular 7r-orbit, 
because the p3/2-shell closure can not be described by simple 7r-orbits, Therefore, in the following discussions, we 
concentrate on the positive-parity orbits of the valence neutrons associated with the molecular cr-orbits and (5'-orbits. 

Fig. [j] shows the density distributions of the single-particle wave functions of the first and second highest neutron 
orbits. In the low-spin cluster states, the positive-parity orbits of the valence neutrons can be well associated with the 
two types of the molecular orbits (cr and 6'). In the state, two valence neutrons with up and down spins occupy 
the cr-like orbits, which have 2 nodes along the longitudinal axis. In the 0^j~ state, which is dominated by the other 
2Tlu configurations, the two neutrons occupy c^'-like orbits. It is very surprising that the developed 6 He+ 6 He cluster 
structure in the O3" state is understood by the new molecular <5'-orbits. It occurs when 2 deformed 6 He clusters are 
attached in parallel. In the 0^ state, all of the 4 valence neutrons are in the negative-parity orbits. Comparing the 
energies of the 0^ state with those of the 0^ and 0^~ states, the 5'-orbit is the highest among the molecular orbits (cr, 
tt and 8'). 

The molecular cr-orbit is one of the reasons for the deformed ground state of 12 Be with the 2huj configurations, which 
is lower than the closed neutron-shell state. Since Be nuclei prefer prolate deformations because of the 2a-cluster 
core, the cr-orbit gains kinetic energy in the developed cluster system. In pioneering studies, |9|]Il|], the importance 
of the cr-orbit in the ground states of n Be and 12 Be were discussed in relation to a vanishing of the magic number. 
Thus, the neutrons in the cr-orbit play an important role in the cluster development of the ground state of 12 Be. On 
the other hand, in case of the 0^" state, we consider that the cluster development is further enhanced due to the effect 
of the orthogonal condition to the 0^ state. 

The next problem is whether or not the molecular orbits appear in the unnatural spin-parity state, which has 
another deformed structure dominated by the 2huj configuration. According to an analysis of single-particle wave 
functions, as shown in Fig. ^ we find that the positive-parity orbits of the 3^ state consist of a cr-like orbit and a 
(5'-like orbit. As mentioned before, the total intrinsic spin of neutrons equals to that in the state (see Table IV). 
This means that the intrinsic spins of two neutrons in the cr-orbit and the J'-orbit are aligned to be one. It is a unique 
character of the 3^~ state to be different from the other 2Uuj bands where the intrinsic spins of two neutrons in the 
same spatial positive-parity orbits couple off to be zero. 

As mentioned above, in the case of the low-spin states in the deformed band, the highest two positive-parity orbits 
of the valence neutrons can be well associated with the cr and <5'-orbits. However the positive-parity orbits in the 
high-spin states can not be classified by the simple a and J'-orbits, but are related to the mixed orbits of the two. For 
example, the highest neutron orbit in the 6+ states indicates the mixing of a and 5', as shown in Fig. ^. As the total 
spin increases from 0^ to 6^ in the K n = 0^ band, the alignment of the neutron intrinsic spins grows accompanying 
mixing the (5'-orbit in the cr-orbit. 

The states in the negative-parity band(if 7r = 1^) are described by 1/kJ configurations, because one of the highest 
2 orbits is dominated by the positive-parity component, while the other one is almost a negative-parity orbit. The 
interesting point is that the 30% mixing of the positive-parity component in the negative-parity orbit gives rise to the 
axial asymmetric shape of the intrinsic state, as shown by the neutron density (Fig. ||). 

In an analysis of 12 Be with the molecular orbits, it should be stressed that all of the states with a 2a core can 
not be necessarily described by pure molecular orbits: tt, a and c^'-orbits. Firstly, the single-particle orbits of the 
valence neutrons are not parity-eigen orbits. For example, although the highest two neutron orbits in the 0^~ state are 
dominated by positive-parity components, they contain 20%^30% mixing of negative-parity contamination, which 
can not be regarded as the 7r-orbits but includes higher-shell configurations. Moreover, mixing of the cr and ^'-orbits 
is found in the high-spin states in the K v = 0\ and K* — Ojj" bands, as can be seen in the highest single-particle 
orbits of the 6+ state (Fig. ||). 



D. Systematics of the 2a clustering in Be isotopes 

When we roughly regard the valence neutron orbits in the Of, lj~, 0j states of 12 Be as the 7r-orbits and the cr-orbits, 
these 12 Be states remind us of analogous states in 9 Be, 10 Be and 11 Be. In terms of the molecular orbits surrounding 
the 2a core, the typical rotational bands in these Be isotopes are classified by the number of occupied cr-orbits. Here, 
we use the notation Tx m u n for those states having m neutrons in the 7r-orbits and n neutrons in the cr-orbits around 
a 2a core. The valence neutron orbits in the lowest natural parity states( 9 Be(l/2"), 10 Be(0f), n Bc(l/2~) and 
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FIG. 7. Relative distance between 2-a clusters in the intrinsic state of the major component described by a single AMD 
wave function. The distance d is evaluated by the relative distance between the mean centers of 2 proton pairs, defined as 
d = y(Re[Zi + Z2 — Z3 — Z4]) 2 [ly/v. The dotted, thick solid and thin solid lines correspond to the distances in the band 
head states with zero(<r°), one(cr 1 ) and two(cr 2 ) neutrons in the a orbits. Namely, the lines (a , a 1 and a 2 ) correspond to 
{ 9 Be(l/27), 10 Be(0+), n Be(l/27), 12 Be(0+)}, { 9 Be(l/2+), 10 Be(17), n Be(l/2+), 12 Be(17)} and { 10 Be(0+), n Be(3/2 2 -), 
12 Be(0+)}, respectively. The open circles indicate the ground states. 

12 Be(0^)) are described as 7r"V , and those in the lowest unnatural parity states ( 9 Bc(l/2+), 10 Be(l~), n Be(l/2+) 
and 12 Be(l~)) are noted as ^'""V 1 , while the 10 Be(0^), n Be(3/27) and 12 Be(fr|) states correspond to n m - 2 a 2 , 
where m is the number of valence neutrons around the 2-a core (see also Refs. |3p, |i~l| , p4| ) . It is important that the 
number of neutrons in the cr-orbits influences the spatial development of the 2a core. The changes in the relative 
distance between the 2 a clusters in the Be isotopes are shown in Fig. (?]. The dotted, thick solid and thin solid lines 
correspond to the distances in the band head states with zero (ir m a°), one {■K m ^ 1 a 1 ) and two (7r m_2 cr 2 ) neutrons 
in the cr-orbit. In each Be isotope, the 2-a distance is smallest in the cr° state and becomes larger and larger with 
increasing number of the occupied cr-orbits, because of the energy gain of the cr-orbits in the developed cluster system. 
Comparing the 2-a distance in each line, the dependence on the neutron number shows that the cluster development 
systematically decreases as the number of the occupied 7r-orbits increases. Especially, a drastic change in the cr° states 
from 9 Be to 10 Be is caused by the P3/2 sub-shell closure effect. Also, in the a 1 and cr 2 lines, there exist gaps between 
the states with a closed p3/2-shell and those with an open p 3 /2-shell. Thus, the closed p3/2-shell plays an important 
role in weakening the 2a development. 

One of the interesting features of neutron-rich Be isotopes is vanishing of the neutron magic number N = 8 in 11 Be 
and 12 Be. In Fig. [j], the ground states are plotted as open circles. The intruder ground states in 11 Be and 12 Be 
correspond to points on the cr 1 and cr 2 lines, respectively. As shown in the Fig. ^, the 2a distances in the ground 
state of 11 Be and 12 Be are about 3 fm. We may conjecture that the magnitude of 3 fm is the natural distance between 
2 a, which is energetically favored in these neutron-rich Be isotopes. Therefore, one of the reasons for the intruder 
ground states can be understood to be a restoration of the natural distance at about 3 fm in the very neutron-rich 
Be isotopes. 

The N = 6 sub-shell closure effects are clearly seen in the energy gaps of 6 MeV as well as in the drastic structure 
changes between the cr° and cr 1 states in 10 Be. The sub-shell effect is also reflected in the structure change between 
the cr 1 and cr 2 states in 11 Be, as shown in the 2a distance. This means that we can regard the vanishing of the magic 
number in neutron-rich Be isotopes as a shift of the neutron magic number from N — 8 to N — 6. 

E. Inter-cluster motion and He decay width 

In order to investigate the inter-cluster motion, we extracted the relative wave functions between clusters in the 
6 He+ 6 He and 8 He+ 4 He channels. We assumed that the intrinsic wave functions of the He clusters are given by the 
+ states of 4 He, 6 He, 8 He in the SU3 limit with the same width parameter as that of 12 Be. The detailed structures 
of the 6 He nucleus, such as the neutron halo, are omitted in the present analysis for simplicity. We estimated the 
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inter-cluster motions between clusters in 12 Be by the reduced width amplitudes which were calculated by 

projecting the 12 Be wave functions to the cluster model space expressed by the superpositions of Brink functions. 
Moreover, the spectroscopic factors (S) and the total cluster probabilities (P c ) are helpful to quantitatively discuss 
the development of He+He clustering. The detailed definitions and the practical calculation of yz,(r), S and P c are 
explained in appendix A. 

Figs. H and ||show the R.W.A. (rj/i,(r)) in the 6 He+ 6 He and 8 He+ 4 He channels extracted from the obtained J = L 
states of 12 Be, respectively. The magnitudes |ryi,(r)| at the channel radius a = 5 fm, the spectroscopic factors S and 
the total cluster probabilities (P c ) are presented in Table ^. 

In each band, even if a low spin state has a large amplitude Iry^ (r) | at the surface region and a large spectroscopic 
factor S, the amplitude becomes smaller in the high-spin states with J > 6. In particular, the amplitudes in the J = 8 
state are very small. We consider that one of the reason for this reduction of the R.W.A. may be the effect of spin- 
alignments in the high spin states. In the 6 He+ 6 He channel, the amplitude at the surface region is most remarkable 
in the Ojj" state, which has the peak amplitude at the surface arround r = 5 f m and a largest spectroscopic factor. It 
indicates the spatially developed 6 He+ 6 He clustering in this state. Compared with that in the 0;j" state, the amplitude 
in the 0^ state shifts inward and has a peak at about r = 4 fm. In the 0^ state, the surface amplitude is smallest 
among three + states, which is consistent with this state has no developed clustering but the dominant neutron 
p-shell closure components. The amplitude |ryx,(r)| at about 4 fm in this 0^" state are caused by the mixing with 
the Of state. Concerning the J ± = 2f ,4^, and 6 J states, both the spectroscopic factors and the surface amplitudes 
\ayL(a)\(a = 5 fm) for the 6 He+ 6 He channel are not small, but those for the 8 He+ 4 He channel are small compared 
with the 6 He+ 6 He channel. On the other hand, with respect to the R.W.A for the 8 He+ 4 He channel, the surface 
amplitudes at r — 4 fm are found to be large in the states in the K v — band, while they are small in other bands 
such as K v = 0^. By analyzing the spectroscopic factors(Table M) and the R.W.A. (Fig. || and ^), we conclude that 
the K~* = 0^ band has both components of 6 He+ 6 He and 8 He+^Ie cluster structures at least in the low spin states. 
It is consistent with the result by the coupled-channel cluster model with 6 He+ 6 He and 8 He+ 4 He wave functions in 
Ref. 0. On the other hand, in the K 71 = Of band, the 6 He+ 6 He clustering appears predominantly. These results 
strongly suggested that the neutrons around the 2a core move about over the whole system in the K v = Of band 
states, while in the K T = Og" band states they move about not over the whole system, but around either of the two 
a clusters. 

We next mention the partial decay widths of the excited states above the threshold energies of the He decays. We 
calculated the theoretical values of the partial decay widths (r4 He and r6jj ) concerning the simple binary decays, 
6 He(0 + )+ 6 He(0 + ) and 8 He(0+)+ 4 He(0+), using the method of reduced width amplitudes: 

r^He = 2P L (a) ■ 74 2 6He (a), (3) 

Pi ( a ) = F*(ka)+Gl(ka)> ( 4 ) 

74 2 6H >) = &^l(«)| 2 , (5) 

(6) 

where k is the wave number of the resonance energy (E r ), k — \J 2[iE r /h 2 , a is the channel radius, and Fl and Gl are 
the regular and irregular Coulomb functions. The resonance energies were evaluated by subtracting the experimental 
threshold energies of the channels from the theoretical excitation energies. In principle, since the present AMD wave 
functions are not sufficient to describe the long tail of the resonance states, we should smoothly connect the tail of 
the irregular Coulomb function (Gi(fcr)) to the original relative wave function ry^ir) for the resonance states at the 
surface region as usually done in the bound state approximation. We found that we can smoothly connect ri/i,(r) 
with Gi(fcr) at the point around r = 5 fm in most of the resonance states in the present case. Therefore we chose the 
channel radius to be a — 5 fm to evaluate the partial decay widths. In Table W% the theoretical values of the reduced 
widths 7 2 He (a), 7 2 H C ( a )' anc ^ partial decay widths are listed. The width of the 0^" state for the 6 He decay is 
broad because of the developed 6 He+ 6 He cluster structure and the lack of a centrifugal barrier. Although the excited 
states seem to be stable for the partial decay width of the He channels, the stability of the resonance states should be 
carefully investigated by taking all of the other possible decay channels, such as the neutron decays and the excited 
He decays, into account. 

VI. SUMMARY 

We studied the structures of the ground and excited states of 12 Be based on the framework of the AMD method. 
This was the first microscopic calculation which systematically reproduced the energy levels of all of the spin-assigned 
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012345678 

Inter-cluster distance r (fm) 

FIG. 8. The inter-cluster motion in the 6 He+ 6 He channel in 12 Be. The reduced width amplitudes(R.W.A), ryL(r) is defined 
in the text and appendix. The thick solid(dashed) lines are H/i(r) of the states in the K n = O3 (K* = 0+) band. 
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FIG. 9. The inter-cluster motion in the 8 He+ 4 He channel in the ground and excited states of 12 Be. The wave function rj/i(r) 
is defined in the text and appendix. 
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TABLE V. Approximated values for the S factor, the total cluster probability of the He-cluster channels in Be. The 
magnitudes, \ayL(a)\, of the reduced width amplitudes at the channel radius a = 5 fm are also shown. 
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TABLE VI. The theoretical values of the partial decay widths and the reduced widths for the 6 He+ 6 He and 8 He+ 4 He decays 
from excited states of 12 Be. The channel radius is chosen to be a = 5 fm. 
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states in 12 Be, except for the 1~ state. One of the present discoveries is the K n = Of band with the 6 He+ 6 He cluster 
structure, which corresponds well to the states recently observed in the He-He break-up reactions. The ground-state 
properties, such as the /3-decay strength and the El transition strength, were reproduced by the present calculations, 
which suggest that the ground state is dominated by a deformed intruder state with the developed cluster structure. 
The theoretical results predicted many low-lying excited states. It was found that most of them have a 2a core, while 
the 1~2 state has a no a-cluster structure. 

We analyzed the structures and the single-particle wave functions of the intrinsic states, and found rotational bands 
with K 7 * = Of, Oj, 1^,03 , 2+, which are dominated by the 2hw, Ofkv, lhu), 2Hlo and 2ftw configurations, respectively. 
The 2hu> and lhu> states have deformed structures with a 2a core. In the positive-parity orbits occupied by the valence 
neutrons, we found the molecular orbits, a-orbits and the new ^'-orbits. The positive-parity molecular orbits play 
an important role in the development of cluster structure. Especially, the reason for the intruder ground state can 
be understood by the energy gain of the er-orbit in the developed cluster state. Although the molecular orbits are 
helpful to understand the structures of the band head states, the model of the 2a and 4 neutrons in the molecular 
orbits is too simple to describe all of the excited states, because such phenomena as the intrinsic spin alignment, the 
dissociation or breaking of 2a core and the mixing of the molecular orbits appear in 12 Be. 

We discussed the inter-cluster motions in the 6 He+ 6 He and 8 He+ 4 He channels, and estimated the partial width 
of these decay channels with the method of reduced width amplitudes. The analyzed results of the partial widths 
strongly suggest that the neutrons of the K* = Of band move around either of two a clusters, while the neutrons of 
the if" = Of band move over the whole system. 

We also discussed the systematics of the cluster development in the neutron-rich Be isotopes according to the 
classification of the states with the number of neutrons in the cr-orbits. Concerning the vanishing of the neutron 
magic number N = 8 in the Be isotopes, we suggested that the reason for the intruder ground states may be 
explained by the restoration of the natural distance between the 2 a clusters. 
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APPENDIX A: INTER CLUSTER WAVE FUNCTIONS 



The calculational methods of the reduced width amplitudes(R.W.A.), spectroscopic factors, and total cluster prob- 
ability concerning the 6 He(0 + )+ 6 He(0 + ) and 8 He(0 + )+ 4 He(0 + ) channels are described in this appendix. 



1. Reduced width amplitudes and cluster model space 

The reduced width amplitudes, are defined as follows: 

y L (a) = g 1 <^^Ybo(O^o(C 1 )0o(C a )|$ i ), (Al) 




ft= t= , \ i h ( A2 ) 



where $£, is the internal wave function of a model wave function and <fio(Ci) and ^0(^2) are the internal wave functions 
of the clusters Ci and C2. The mass numbers of the system and the clusters Ci and C2 are A, A\ and A2, respectively. 

Here, we briefly review the R.W.A. in the system of two clusters written in the form of RGM(resonating group 
method) or GCM (generator coordinate method) wave functions. The detailed calculational methods of the RGM and 
GCM kernels are described, for example, in Rcf. |53|. We assume the same width parameter (y) for clusters Ci and 
C2 described by harmonic oscillator shell model wave functions. When $l is an RGM- type wave function expressed 
as, 

*l = q2A{ X L(r)Y L0 (f)MGi)MG2)}, (A3) 
92 = -J L_ _, (A4) 

the R.W.A. , 2/i(a), can be calculated based on the knowledge of the RGM norm kernel as follows. By expanding the 
relative motion xl{t) with the radial harmonic oscillator (H.O.) wave functions R n i{r, v 1 ) with the width parameter 



^*/as 



Xl(t) = y £ d e nL R nL {r,v'), (A5) 
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we obtain the R.W.A. as, 

VL{a) = ^ e nL^nLRnL{a, v'), (A6) 



where /j, n L are the eigen values of the RGM norm kernel. If the RGM-type wave function A3 is normalized to unity, 
the following equation is satisfied: 

Y,e 2 nL ^nL = l- (A7) 

n 

In the case of general A-nncleon wave functions, $l can be separated into the cluster part and the non-cluster part, 

*l = A{ X L(r)Y L0 (f)MGi)MG2)} + (A8) 

where is the residual part after projection of the cluster model space, and satisfies (i / Lo(^)0o(Ci)^o(C2)| < i ) f ) = 0. 
In this case, the R.WA. of <I>l can be obtained with equation |A6| by using an expansion of Xi( r ) with the H.O. 
functions shown in EqA5. The normal spectroscopic factors S are calculated as 

S = [ y L (r) 2 r 2 dr = J2^nLel L - (A9) 

It is useful to consider the total probability of the clustering component (total cluster probability) in $l defined by 
the squared overlap with the cluster model space calculated as, 

P c = J2^Le 2 nL - (A10) 



It is easily found that P c = 1 for the normalized RGM-type wave functions from equation A7. 
The R.W.A. of the wave function $i( 12 Be) in the AMD framework is 

y L (a) = 9l (fc^y LO (f)0 o (C 1 )0o(C 2 )|$ L ( 12 Be)), (All) 

where $i( 12 Be) is a normalized spin-parity projected AMD wave function, where the center of mass motion is 
extracted, and 0o(C„) (C„ = 4 He, 6 He, 8 He) are the intrinsic wave functions of the SU^Ahait + states of the He 
clusters. In order to project to the cluster model space from the 12 Be w ave functions (<I>i( 12 Be)) and calculate 
the radial wave functions (xi( r )) between the clusters, as mentioned in EqA8, we use the orthonormal sets of the 
GCM wave functions. Concerning the fc-th cluster wave functions ^> c {hk) with the inter-cluster distances we 
adopt Brink-type functions, where two clusters(Ci = Al He and C2 = 2 He) written in the SU3-lim.it representation 
are located at the points (0,0, A ^ A &fc) and (0,0, — A ^ A bk), respectively, 

* c (b fc ) =^{*(Ci > ^b fc )*(C 2 ,-^b fc )}, b fc = (0,0,6 fc ). (A12) 

The width parameters of the He clusters are chosen to be the same as that of the 12 Be wave functions for simplicity. 
Then we can rewrite \E' c (bfe) as 

* c (b fe H^ (X G ).$ c (b fc ), (A13) 
$ c (b fc ) =g 2 ^{r(r,b fe X)<MWo(C 2 )}, (A14) 
T(r, b fc , V') = (%:)3/4 e -'(r-b fc ) 2 (A15) 

co(x G )^(^) 3/ V^ x2 o, x G EiE i=l/i) (Ai6) 

In each L, we make a set of the orthonormal basis <tk,L by linear combinations of the total-spin-projected GCM 
wave functions $>k,L'- 

$ fc , L = g fcL P L $ c (b fc ), (A17) 

*fc,£ = E fc '^*fc',£. (A18) 
($k,L\$k>,L) = 8 k u, (A19) 

(A20) 
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where quh are the normalization factors qkL = l/(^o6$ c (bOl^oo^ c (' :, fc)) 1/ ' 2 - By using the partial wave expansion of 
the functions r(r, bfe,z/), the radial functions Xl ( r ) m * ne ^"^ n wave functions <&k,L can be expressed as follows: 

$fe,L = ! ^{xi fc) (r)r L o(r)^o(C 1 ) ( /)o(C 2 )}, (A21) 



r 



X^M-^V^r^r^,^'), (A22) 

, \ 3/4 



(r,& fe ,z/)= (^) ^i L {2y'rb k )e^'^+ b D, (A23) 

where il is the modified spherical Bessel function. We assume that the projection operator P£ onto the cluster model 
space can be written by the orthonormal basis 



pz = 52\$k,L)($ k , L \= E i*^>4M^*<<^i. (A24) 



In this case, the radial function xl{t) in Eq.|A§| for <f>L( 12 Be) is written as 

,(*') 



Xdr) =Y^x { l \r){J24 L J4 L J*(^,L\^ L C 2 Be)}}. (A25) 



2. Practical calculation of fi n L 

In the present calculations, we chose bk = 1, 2, • • • , 9 fm for L = 0, 2 (L = 0, 2, 4, 6), bk = 2, 3 • • • , 9 fm for L = 4,6, 
and &fc = 3, • • • , 9 fm for L = 8 (L = 8) for the 8 He+ 4 He ( 6 He+ 6 He) channel. In the practical calculations, we 
approximate fi n L with /i„x calculated as follows. It is known that the eigen values ^ n L equal to zero for the Pauli 
forbidden states with N = 2n + L < N rn i n {N m i n is the minimum allowed number) and /i„l ~ 1 with enough large 
numbers, N = 2n + L. In order to obtain approximated values we truncate the quanta n with a finite number 
n = n mm , ■ ■■ , n max (N rnin = 2n min + L, N max = 2n max + L) by assuming fi nL = 1 for N = 2n + L > N max . With 
the use of the expansion of x^l ( r ) m ^ ne normalized cluster wave functions <f>k,L by RnL,(r, v 1 ) 



(A26) 



we obtain /i n £ by following equations: 



E(fe)2 - -. 
e nL finL = 1. 



(A27) 



The approximated values p, n L for simple systems, a + a and a+ 16 are compared with the exact eigen values /j, n L in 
Table IvIB. 



In the present calculation of 12 Be, N max was chosen to be 14. By using /2 n £ determined by Eq.A27, we define the 
approximate R.W.A., yx,( a), t he spectr osco pic facors, S, and the total cluster probability, P c for $i( li! Be) by analogy 
with the relations in Eqs. A6, AG and AlCj : 



VL{a) = J2 n e n Lfl>nLRnL{a,v'), 



(A28) 
(A29) 
(A30) 



where the coefficients e„L were calculated by e„L = J r 2 XL(r)R n Ldr from Xi( r ) given in Eq. A25 
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TABLE VII. The eigen values ^N=2n+L and approximated values /1„l(L — 0) of the RGM norm kernel for a + a and a+ 16 
systems. The values fijsf are taken from Ref. [H. The width parameters for a and 16 are assumed to be same. 
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